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Abstract 

f2mma program can be used for translation a program written in some subset of the 
FORTRAN language into Mathematica programming language. This subset have 
been enough to translate GAPP (Global Analysis of Particle Properties) program 
into Mathematica language automatically. A table with Standard Model observables 
calculated with GAPP (Mathematica) is presented. 
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Nature of problem: 

f2mma program allows to translate programs from some subset of FORTRAN pro- 
gramming language to Mathematica system's language (where mathematical objects 
can be analyzed more efficiently in symbolic form). The general aim of f2mma is 
translation of GAPP (Global Analysis of Particle Properties) package [1]. 

Solution method: 

Perl::RecDescent module was used with corresponding FORTRAN grammar and 
transformation rules. 

Unusual features: 

f2mma can generate additional new lines in Mathematica output programs. 
Additional comments: 

(1) f2mma use only some subset of FORTRAN language grammar (for example, 
f2mma goes not process 'goto' operator); (2) f2mma does not make syntax check of 
input FORTRAN program; (3) translation process can be slow for large FORTRAN 
programs. 

References: 

[1] J. Erler, hep-ph/0005084 
[DOWNLOAD 'f2mma' (1.2 Mb)| 
LONG WRITE-UP 

1 Introduction 

In programs written in FORTRAN programming language scientific models 
of specific field have been usually embedded into codes. It's difficult and im- 
practical to support or to upgrade them especially by non-programmers. So 
the only aim of the codes is often to carry out numerical calculations. 

From other side, it's known that with help of computer algebra system such as 
Mathematica one can represent the scientific model as an object which can be 
simply analyzed and improved. Mathematica provides important capabilities 
such as: 

• High-precision calculations; 

• Special functions; 

• Equations solvers; 

• Matrices and vectors manipulations; 

• Sophisticated programming language; 
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• 2D and 3D graphics. 

So, there is an idea to translate FORTRAN codes to Mathematica system's 
language and then carry out some calculation in symbolic form. f2mma pro- 
gram can be used as such translator. 

For others motivations see [5], for example. 

A design of f2mma program had not aim to provide tool for translation every 
FORTRAN program (according to any FORTRAN language standard) into 
Mathematica language but actually to translate GAPP 1 [1] package into CAS 
language (it's important for us since GAPP contains recent scientific models 
for electroweak interactions of the particles). This purpose was achieved. In 
the section 3 we present a table with Standard Model observables calculated 
with GAPP (Mathematica). 

f2mma is a rather simple program and can be easily extended by user. Special 
effort was made to document the codes. 



2 An example 

Let's translate a test program (which should print out several 'OK' messages 
and expressions equaling to zero) in order to see the basic functionalities of 
f2mma program: 

TEST RUN INPUT 

c Test subroutine 

subroutine subl(x) 

c ... types of scalar variables are ignored: 

integer x,y 

common /t/ r 

real*8 r 

x=x+l 

return 

end 

c Test function fl 

integer function fun(x,y) 
integer x,y 

integer ml (2), m2(2:4), m3(x:y) 
c ... test 'data' 

1 GAPP - Global Analysis of Particle Properties 
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data ml /ll, 22/ 

print *,'ml(l) - 11 = ' , ml(l)-ll 

print *,'ml(2) - 22 = ', ml (2) -22 

m3(x)=x+y 

fun = m3(x) 

return 

end 

c 

c ... '_' symbol is replaced by > TTT' string: 
integer x, y, z, t, ar_int (1 : 5,2 : 5) , n, si, s2, fun 
common /t/ r 
real*8 rl,r2 

double precision dpi, dp2, dp3, dp4, prec 
complex*16 q 
logical flag 

print *, '6 OK messages should appear:' 

x=l 

rl=l. 

r2 = 1 . 3 

c ... decimal fraction converts to expressions: 

dpi = 1/1. d3 

ar_int(l,2) = 1 

dp2= 1 . 123 456 d-7 

flag = .true. 

prec =l.d-100 

c Test translation of mathematical expressions: 

dp3=- (l+2**x/2) /dsqrt (y**2+3 . dO) +1 
if (dp3 - prec.le. 0) print *,'0K-1' 
z=-x**2 
y=(-x)**3 

if (z-y.eq.O) print *,'0K-2' 
t = z-y 

c Complex numbers : 

q=(0,l)**2+(l,0)**2 
dp4=cdabs (q) 

if (dp4 - prec .le. 0) print *,'0K-3 ; 

c Test complex expressions with 'if: 

if (.true.) then 
print *, '0K-4' 
else 

if (.true.) then 
t=l 

else 

t=2 

endif 
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endif 

if (t.eq.O .and. f lag. eqv. . true . ) then 

print *, 'OK-5' 

else if (.true.) then 

t=4 

else 

t=5 

endif 

n=9 

sl=0 

c Test 'do' : 

do 10, x=0,n 
10 sl=sl+x 
s2=0 

do 20 x=n,0,-l 
s2=s2+x 

20 continue 

if (sl.eq. (n*(n+l)/2) . and. si . eq. s2) print *, 'OK-6' 
x=0 

y=o 

r=0 . eO 

print *, 'Four on right-hand side should appear:' 
call subl(n) 

print *, 'n - 10 = ' , n-10 
n=fun(l,2) 

print *, 'n - 3 = ' , n-3 

end 

TEST RUN OUTPUT 



(* c Test subroutine *) 

Set Attributes [subl, HoldAll] ; 
subl [x_] :=Module [{ y }, 

(* c ... types of scalar variables are ignored: *) 

(* #common: r *) 

x=x+l; 

Return [subl] ; 

]; 

(* c Test function fl *) 

Set Attributes [fun, HoldAll] ; 
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fun[x_,y_] :=Module [{ fun,m3,ml,m2 }, 



F2MmaDimensions[ml]={{l, 2}}; 
F2MmaDimensions[m2]={{2, 4}}; 
F2MmaDimensions [m3] ={{x, y}}; 

(* c ... test 'data' *) 

(* # Inserting to ml *) 
Module [{F2MmaTempl, F2MmaTemp2, j, jO, k}, 
F2MmaTempl={ 11,22 }; 
F2MmaTemp2=F2MmaDimensions [ml] ; 

j=Array [jO, Length [F2MmaDimensions [ml] ] ] ; j=Transpose [F2MmaTemp2] [[1]] ; 

Do [k=l ; ml [Sequence@@j] =F2MmaTempl [ [i] ] ; 

j[[k]]++; 

While [k<Length [F2MmaTemp2] && j [ [k] ] >F2MmaDimensions [ml] [ [k , 2] ] , 
j [[k]]=F2MmaDimensions[ml] [[k,l]] ; j [[k+l]]++; k++] ; 
, {i , 1 , Length [F2MmaTempl] }] ; ] ; 

Print ["ml(l) - 11 = " ,ml [1] -11] ; 

Print ["ml (2) - 22 = " ,ml [2] -22] ; 

m3 [x] =x+y ; 
f un=m3 [x] ; 
Return [fun] ; 

]; 

(* c *) 

(* c ... '_' symbol is replaced by 'TTT' string: *) 
F2MmaDimensions[arTTTint]={{l, 5}, {2, 5}}; 
(* #common: r *) 

Print["6 OK messages should appear:"]; 

x=l; 
rl=l; 

r2=(13*10~(-D) ; 

(* c ... decimal fraction converts to expressions: *) 

dpl=l/(l*10~(3)); 
arTTTint[l,2]=l; 
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dp2=(1123456*l(r (-6)*1(T (-7)) ; 
f lag=True; 
prec=(l*10~(-100)) ; 

(* c Test translation of mathematical expressions: *) 

dp3=-(l+2~x/2)/Sqrt [y"2+ (3*10" (0) )] +1 ; 
If [ dp3-prec<=0 
, Print ["0K-1"] ; 

]; 

z=-x~2; 

y=(-x)~3; 

If [ z-y===0 

, Print ["OK-2"] ; 

]; 

t=z-y; 

(* c Complex numbers: *) 

q=(0+I*(l))~2+(l+I*(0))-2; 

dp4=Abs [q] ; 

If [ dp4-prec<=0 

, Print ["OK-3"] ; 

]; 

(* c Test complex expressions with 'if: *) 

If [ True, Print ["0K-4"] ; 

, If [ True, t=l; 
, t=2; ]; 

]; 

If[ t===0 && flag===True, Print ["OK-5"] ; 

, If [ True, t=4; 
, t=5; ]; 

]; 

n=9; 
sl=0; 

(* c Test 'do' : *) 

For[x=0,x<=n,x+=l,sl=sl+x; ]; 
s2=0; 
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For[x=n,x>=0,x+=-l,s2=s2+x; ]; 
If [ sl===(n*(n+l)/2) && sl===s2 
, Print ["OK-6"] ; 

]; 

x=0; 
7=0; 

r=(0*10~(0)) ; 

Print ["Four on right-hand side should appear:"]; 
subl [n] ; 

Print ["n - 10 = " ,n-10] ; 
n=fun[l,2] ; 

Print ["n - 3 = " ,n-3] ; 



3 Table with Standard Model observables 



In order to check general correctness of f2mma translation we present ta- 
ble with SM observables (that is values of physical expressions arising from 
Standard Model theory of particle physics). They were calculated with fallow- 
ings values of the fitted parameters: M z = 91.187 GeV, m t = 166.96 GeV, 
m b = 4.21 GeV, m c = 1.29 GeV, a s (0) = 0.121, m H = exp(4.7126) GeV, 
T = U = S = B = 0, Z = 3, M z > = 1000 GeV, sin9 = 0, X g = 1. Total 
X 2 (with parameters defined above) is about 49.6 for 37 effective degrees of 
freedom (it's quite comparable with G APP (FORTRAN) result 49.4 at the 
minimum of x 2 )- 

For explanation of the observables abbreviations please see [1], [2]. 
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Table 1 



Table of the results obtained with GAPP (Mathematics). 'SM value' - calculated 
value, 'Exp.' - experimental value of observable, 'Exp. una' - experimental uncer- 
taincy, 'Pull' = ('SM value' - 'Exp.')/'Exp. unc' 



ID 


Observable 


SM value 


Exp. 


Exp. unc. 


Pull 


1 


M z 


91.187 


91.1876 


0.0021 


-0.285714 


2 


Tz 


2497.06 


2495.2 


2.3 


0.807306 


3 


Chad 


41.4737 


41.541 


0.037 


-1.81758 


4 


Re 


20.7482 


20.804 


0.05 


-1.11679 


5 


Rfj, 


20.7483 


20.785 


0.033 


-1.11132 


6 


R T 


20.7936 


20.764 


0.045 


0.65709 


7 


A FB (e) 


0.0162736 


0.0145 


0.0025 


0.709421 


8 


A FB (mu) 


0.0162736 


0.0169 


0.0013 


-0.481883 


9 


A FB {tau) 


0.0162736 


0.0188 


0.0017 


-1.48615 


10 


P(tau) 


0.147303 


0.1439 


0.0043 


0.791313 


11 


P FB {tau) 


0.147303 


0.1498 


0.0049 


-0.509664 


12 


sin 2 (9 e e ff)(Q FB ) 


0.0422812 


0.0403 


0.0026 


0.762011 


13 


A FB (s) (DELPHI + OPAL) 


0.10337 


0.0976 


0.0114 


0.506097 


14 


Rd,s/ (Rd + Ru + R s ) 


0.359181 


0.371 


0.023 


-0.513853 


15 


Rb 


0.215635 


0.21638 


0.00066 


-1.12942 


16 


Rc 


0.172329 


0.172 


0.003 


0.109588 


17 


A FB (b) 


0.103262 


0.0997 


0.0016 


2.22614 


18 


A FB (c) 


0.0737844 


0.0706 


0.0035 


0.909841 


19 


A LR FB (b) 


0.934691 


0.925 


0.02 


0.48454 


20 


A LR FB (c) 


0.667872 


0.67 


0.026 


-0.0818604 


21 


Aiji(hadrons) 


0.147303 


0.15138 


0.00216 


-1.88766 


22 


ALFiileptons) 


0.147303 


0.1544 


0.006 


-1.18289 


23 


A LR FB {mu) 


0.147303 


0.142 


0.015 


0.35351 


24 


A LR FB (tau) 


0.147303 


0.136 


0.015 


0.75351 


25 


A e (Q LR ) 


0.147303 


0.162 


0.043 


-0.341799 


26 


A LR FB (s) 


0.935665 


0.895 


0.091 


0.446873 


27 


M W {LEP) 


80.391 


80.412 


0.042 


-0.498945 


29 


M\y(Tevatron) 


80.391 


80.45 


0.058 


-1.01648 


30 


Ty/{Tevatron) 


2.09364 


2.103 


0.106 


-0.0883258 
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Table 2 



Continue of the table with results obtained with GAPP(Mathematica). 



31 

o 1 


rnAr)n]p\('\ CDF 1\ 

1 1 if \jJyJl t- J I ± ^ i_s 1 l j 


1 77 007 

J. ( 1 .Uu 1 


1 7fi 1 


7 3fi 

i .OU 


19399R 


39 


mJnnlp\C) CDF T\ 


1 77 007 


1fi7 4 

J-U 1 . t: 


1 1 39 


D 843499 


oo 




1 77 DD7 


1 8fi 


1 1 ^1 


D 781 981 

U. ( OliOl 




m±(nnlf>\('\ DO T\ 


1 77 DD7 

J- f 1 .UU l 


18D 1 


^ 14 


D ^1791 98 




m + (nn1p\( c ) DO T) 


1 77 007 


1fi8 4 


1 9 84 


R703R3 

U.U ( uouo 




1 1 if [jJUi c M X , w is r J- J- J 


1 77 007 

J. i ( .UU 1 


1 77 ^ 


±o. xu 


-u.uo I '-too I 


37 


mAvole)(2 CDF II) 


177.007 


175. 


18.95 


105934 


43 


777 ~ { 711 ~ ) 


1.48297 


1.46484 


0.17046 


106346 


44 


771 L I TH L 1 


1.96517 


1.97361 


162013 


-0 0521324 


45 


Aa^ /(l 8GeV) 


0058189 

KJ • KJ \J t^f 1~ KJ is 


005768 


0.0001 


0.509 


46 


f — 2 — alir)/2 


4509.3 


4511.07 


0.8 


-2 2071 


47 


-tt, T 


292.3 


290.87 


0.52 


2 75057 


48 


rn 2 ( NuTpV 2002) 


303984 


30005 


001 37 

U .UU IJ 1 


2 871 1 7 


4Q 


ar 2 ( NuTpV 2002) 


0300695 

U.UOUUUi7U 


03076 

U. UOU 1 u 


001 1 

U .UU J. J. 


-0 627763 

V / . U Z. 1 1 U(J 


50 


karma(CCFR 1997) 


583381 


0.582 


0.0041 


336865 


51 


7? (CHARM 1984) 


309273 


0.3021 


0.0041 


1.74949 


52 


/? (CDHS 1984) 


309273 


3096 


0.0043 


-0 0760666 

U.U ( uuuuu 


53 


P-fr* ff 4 P /If 1 984) 


38R309 


403 


01 fi 

U .U J-U 


-1 0439 


54 


£>-(CDHS 1984) 


38R309 


384 

U. OOI 


01 8 

u .u ±o 


1 989fifi 

U. l^O^UU 


UU 


p-tanHq iq7q) 


U.OOl / Z<U 


^fi^ 
u.ouu 


01 

U.U J-U 


1 04^4 


UU 


m/ (^ e ) (CH ARM TT) 

u y \ 1 1 J±l Ll vl ll) 


-0 0392145 


-0 04 


015 

U .U J- u 


0523635 


57 


g A {v ' e) {CHARM II) 


-0.506548 


-0.507 


0.014 


0.0323025 


58 


g ee (SLAC £158) 


0.0441208 


0.059 


0.0122 


-1.2196 


60 




-73.2353 


-72.84 


0.46 


-0.859329 


61 




-116.879 


-116.4 


3.64 


-0.131507 


62 


In 73(6 -► sj)/B(b -» cei/) 


-5.73228 


-5.69 


0.17 


-0.248685 


63 


A FB {e){CDF II) 


0.231486 


0.2238 


0.005 


1.53724 
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